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£ I Abstract 

fN| ' Simple, self-similar, elliptic solutions of non-relativistic fireball hydro- 

dynamics are presented, generalizing earlier results for spherically sym- 
' metric fireballs with Hubble flows and homogeneous temperature profiles. 

The transition from one dimensional to three dimensional expansions is 
1 investigated in an efficient manner. 

Introduction. Recently, a lot of experimental and theoretical efforts have 
i| gone into the exploration of hydrodynamical behavior of strongly interacting 

hadronic matter in non-relativistic as well as in relativistic heavy ion collisions, 
see e.g. B - pi. Due to the non- linear nature of hydrodynamics, exact hy- 
dro solutions are rarely found. Those events, sometimes, even stimulate an 
essential progress in physics. One of the most impressive historical example 
is Landau's one-dimensional analytical solution (1953) for relativistic hydrody- 
namics [| that gave rise to a new (hydrodynamical) approach in high energy 
physics. The boost-invariant Bjorken solution [Q, found more than 20 years 
later, is frequently utilized as the basis for estimations of initial energy densities 
in ultra-relativistic nucleus-nucleus collisions. 

The obvious success of hydrodynamic approach to high energy nuclear colli- 
sions raise interest in an analogous description of non-relativistic collisions, too. 
The first exact non-relativistic hydrodynamic solution describing expanding fire- 
balls was found in 1979 ||. It has been generalized for fireballs with Gaussian 
density and homogeneous temperature profiles as well as for fireballs with 
arbitrary initial temperature profiles and corresponding, non-Gaussian den- 
sity profiles. All of these solutions have spherical symmetry and a Hubble-type 
linear radial flow. However, a non-central collision has none of the mentioned 
symmetries. The purpose of this Letter is to present and analyze hydro solu- 
tions for such cases. The results presented in this paper may be utilized to 
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access the time-evolution of the hydrodynamically behaving, strongly interact- 
ing matter as probed by non-central non-relativistic heavy ion collisions ^ |9| . 
As the hydro equations have no intrinsic scale, the results are rather general in 
nature and can be applied to any physical phenomena where the non-relativistic 
hydrodynamical description is valid. 

Generalization of spherical solutions to elliptic flows. Consider an ideal 
fluid, where viscosity and heat conductivity are negligible, described by the local 
conservation of matter (continuity equation), the local conservation of energy 
(energy equation) and the local conservation of momentum (Euler equation): 

dn , 

•V(vn) = 0, (1) 



dt 

de 

Ot 




V(ve) = -PVv, (2) 



mn(-+vV)v = -VP, (3) 

where m is the mass of a single particle and n — n(t, r), v = v(t, r), P = P(t, r) 
and e = s(t,r) are the local density of the particle number, the velocity, the 
pressure and the energy density fields, respectively. To complete the set of 
equations we need to fix the equations of state. For reasons of simplicity we have 
chosen the equations of state to be those of an ideal, structureless Boltzmann 
gas: 

e(t,r) = ^P(i,r), (4) 

P(t,r) = n(t,r)T(t,r). (5) 

The solutions that are presented in the subsequent parts cannot be trivially 
generalized to any arbitrary equations of state. Nevertheless, they provide a 
transparent insight into the collective physical processes in a non-central heavy 
ion collision. 

In the following, we utilize the above ideal gas equation of state and rewrite 
the hydrodynamical equations in terms of the three functions n,v and T. 

In ref. || Q], special classes of exact analytic solutions of fireball hydrody- 
namics were found assuming spherical symmetry and self-similar Hubble flows. 
In ref. || a homogeneous temperature profile was assumed, while the general 
solution for arbitrary, inhomogeneous initial temperature profiles was found in 
ref. 0. In these articles, the concept of self-similarity meant that there is a 
typical length-scale of the expanding system R = R(t) so that all space-time 
functions in the hydro equations are of the form F = G(t)H(s) where s — r 2 / R 2 
is the so-called (dimensionless) scaling variable. 

Let us go beyond spherical symmetry and consider three typical lengths 
of the expanding system: X, Y and Z, all functions of time only. Let us ro- 
tate our frame of reference to the major axis of the ellipsoidal expansion, and 
leave to future applications to relate these major axes to the laboratory frame. 
Consequently, let us introduce three scaling variables x = r 2 /X 2 , y = r 2 /Y 2 
and z — r 2 /Z 2 and assume that all space-time functions are of the form of 
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F = G{t)H{x)K{y)L{z). Using this ansatz we find that the continuity equation 
is satisfied regardless of the density profile if the velocity field is a Hubble-flow 
field in each principal direction: 

v x(t,r) = -TT7-r x , 



v v (t, r) 



X(t) 
Y(t), 



Y(t) 

v z (t,r) = || r ,. (6) 

Although our sole assumption concerning the temperature was the ansatz form 
already mentioned, we found that the Euler equation requires the temperature 
to be homogeneous, independent of the coordinate variables: T = T(t). The 
energy equation is only satisfied if 



T (t)=T [-?-) , (7) 



2/3 

v(t)J 

where V(t) = X{t)Y(t)Z{t) is the typical volume of the expanding system, 
while Vo = V(to) and T = T(to) are the initial temperature and volume. The 
homogeneity of the temperature and the Euler equation implied that the density 
profile is a product of three Gaussians, with different, time dependent radius 
parameters: 

n(t ' r) = no w cxp (-^kf - -ww " ' (8) 

where no = n(0, 0) can be expressed by the total number of particles (N) as 
1 N 

The time evolution of the scales are determined - through the Euler equation - 
by the equations 

T / V \ 2 / 3 

XX = YY = ZZ = — [ — ) . (10) 
m \ V J 

This system of non-linear, second-order ordinary differential equations has a 
unique solution for the scale-functions if the initial parameters Xq, Yq, Zq and 
Xq, Yq, Zq are given. Although this solution has not yet been found in an 
explicite, analytic form, some of its properties are determined in the subsequent 
parts. 

Properties of the elliptic solutions. Global conservation laws reflect, in 
general, boundary conditions for solutions, or their behavior at asymptotically 
large distances. Because of reflection symmetry of densities and velocities, the 
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conservation of momentum, P = 0, is satisfied automatically and gives no non- 
trivial first integral. On the other hand, the asymptotically fast decreasing of 
the densities give us the possibility of using total energy conservation as 

d f ,<> , nrav 2 . ,„ „ . 

m ] + = * (ID 

to find the first integral of the system of equations (0) - (||)- Substitute (ft), (|sj) 
into (O), we get: 



. 2 . 2 . 2 T / V \ 2/3 

X + y + Z +3— = A = const. (12) 



m \ y 



Using (10) one can rewrite (|12|) in the form 



l^(X 2 +Y 2 + Z 2 )=A, (13) 
and find finally 

R 2 (t) = X 2 {t) + Y 2 {t) + Z 2 (t) = A(t - t Q ) 2 + B{t - t ) + C, (14) 

where 

.2 .2 .2 Tn 

A = X Q + Y Q + Zo+3— , 

m 

s = 2pc x +y y +^0 zo), 

C = X 2 +Y Q 2 + Z 2 . (15) 

The simple equation Jl^) express the general property of elliptic hydrody- 
namic flows. The value of radius- vector evolves in time similar to a "particle" 
with coordinates (X, Y, Z) that moves in a non-central, repulsive potential ac- 
cording to eqs. ( J10| ) . In particular, one may introduce the canonical coordi- 
nates (X,Y,Z) and the canonical momenta as (P X) P v , P z ) = m(X ,Y , Z) and 
the Hamiltonian H as a rewritten form of eq. (|12|): 

The Hamiltonian equations of motion can be written in terms of Poisson brackets 
as X = {X, H}, ... , P x = {P X ,H}, ... . The Lagrangian form of these equations 
is given by eqs. (|l0|). Due to the repulsive nature of the potential, the coordinates 
{X, Y, Z) diverge to infinity for large times. As the potential vanishes for large 
values of the coordinates, the canonical momenta tend to constant values for 
asymptotically large times. 

Eq. ( |l2| ) expresses the conservation of whole energy (kinetic and potential) 
of the "particle" , corresponding to H(X, P x , ...) = E = const. The resulting 
eq. (fi~3) has also great importance for the analysis of approximate analytical 
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m / ■ 2 
~2 \ Xas 



solutions. It is worth mentioning another interesting relation that one can get 
from ( |l2] ) for asymptotic times t as , when V(t as ) ^> Vq: 

■2 .2 \ m ( ■ 2 -2 .2\3 , . 

Y as + Z as ) = -[X + Y + Z Q ) + -T . (17) 

This relation expresses the equality of the initial flow and internal energy with 
the asymptotic energy which is present in the form of flow. 

Although eqs. ( |io| ) are easy to handle with presently available numerical 
packages, we note that a further simplification of these equations to a non- 
linear first order differential equation of one variable is possible, if an additional 
cylindrical symmetry is assumed, corresponding to X(t) — Y(t). One may 
introduce the angular variable tj> = arccos(Z/i?) so that 

X(t) = Y{t) = -=R(t)sui(f>(t) (18) 
v2 

Z{t) = R(t) cos <f>(t). (19) 
The time evolution of 4>(t) is determined by the following first order equation: 
1 



R 2 (t) 



E- R 2 {t )-- To{x ° Zo)2/3 



(20) 



m 2 m R 2 (t) (sin <^>) 4 / 3 (cos 

where R(t) is given explicitly by eq. (|l4|). 

Figs. |l]^J indicate the results of numerical solutions of eqs. ([io|). Using 
Landau-type initial conditions, one confirms that even the full three-dimensional 
solution results in a small amount of transverse flow generation, while for more 
general initial conditions, significant amount of transverse flow can be generated. 
Transverse flow is stronger if the initial conditions are closer to spherical sym- 
metry or, if the fraction of the initial thermal energy is increased as compared 
to the initial kinetic energy. For more details, see the figure captions. 

In the last part an approximate, analytic solution is presented that corre- 
sponds to Landau-like, one dimensional expansions. 

Approximate one- dimensional solutions. Hydrodynamical evolution, which 
is described by Eqs. (|l|) - (g), starts from some initial conditions. Consider 
Landau- type initial conditions (longitudinally compressed "ellipsoid"), corre- 
sponding to the real situation in non-relativistic heavy ion collisions: 

Xo=Y = 0, Z « Xo, Zo « Y , (21) 

and in general case Zo^ 0. The last reflects the situation when a system is 
(locally) thermalized before, after or at the moment of full nuclear stopping and 
transverse expansion starts to develop only after the local thermalization. Then 
in some time interval to < t < t the hydro evolution is quasi - one - dimensional: 

f>*i,^;~ (22) 

and the equation of motion for Z(t) takes the following form 
T / Z ^ 2 ' 



Z Z = -± ( ~f ) ■ (23) 
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Figure 1: Top panel shows the time dependence of the major axis of the 
ellipsoidal (Gaussian) density profile, eq. (p|), while the bottom panel indicates 



the velocity of expansion at the rms radii, 



X(t), r y = Y{t) and r z = Z{t), 



for three dimensional, self-similar elliptic flows. Solid lines stand for X{t) and 
v x (t) 1 dashed lines for Y(t) and v y (t), short-dashed lines for Z(t) and v z (t). The 
initial conditions are: Xq — 4 fm, Yq = 3 fm, Zq = 1/4 fm, the initial velocities 
are all vanishing, To/m = 0.1 . Due to the Landau- like initial condition (strong 
initial compression in the z direction) , the flow is almost one dimensional, only 
small amount of transverse flow is generated in this case. 
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Figure 2: Same as Fig. [lj but for a more spherical initial profile with an 
initial inwards flow in the z direction. The initial conditions are: Xq = 4 fm, 
Yq = 3 fm, Zq = 3 fm, the initial velocities are Xq — 0, Yq = 0, Zq = —0.5, 
while To/m = 0.1 . Due to the deviation from the Landau-like initial condition, 
the final flow is almost spherically symmetric, three dimensional, and a large 
amount of transverse flow is generated in this case. 
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Figure 3: Same as Fig. [|, but for To/m = 0.01. This corresponds to increasing 
the initial kinetic energy as compared to the internal (thermalized) energy, and 
in this case the flow becomes approximately one dimensional again. The initial 
conditions are: Xo = 4 fm, Yq = 3 feu, Zq = 3 fm, the initial velocities are 
X = 0,Y = 0, Z Q = -0.5, while T a /m = 0.01 . 
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Figure 4: A more general initial condition may result even in a dominantly 
two-dimensional expansion. The initial conditions are: Xq = 4 fm, Yq = 3 fm, 
Z(0) = 3 fm, the initial velocities are Xq = —0.2, Yq = 0.1, Zq = —0.5, while 
T /m = 0.05 . 
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This equation has an exact analytic solution, 



Z 2 (t)-(Z 2/3 + ( fl+ +a_) 2 ) 3 



(24) 



where 



Z(t) 
Z 



-Z{t)±[Zl + -Z\t) 



W{t-t ), u 2 = -[Z 



1/3 



.To 



T 



3/2 



m 

z o Z fu 2 + To 
m 



3u 4 



(25) 



Here to is the turning point for Z 2 (t) if Zo< 0. Using ( |l4| ) we obtain that the 
conditions for validity of the solution (|2^), (|24|) are satisfied within some time 
interval to < t < t if 



A(t-t ) 2 + B(t-t Q ) + Z 2 - Z 2 (t) 



X 2 o+Yo 2 



< 1. 



(26) 



The hydrodynamical evolution described by the equations (|]J) - (|^) cannot 
be continued infinitely in time, because the general criteria of applicability of 
hydrodynamical description are violated: the mean free path I oc I /(an) has 
to be (much) smaller than the typical length scales of the system, for example 
the effective geometrical sizes or hydrodynamic lengths. Due to the hydrody- 
namical expansion the density ^ will decrease with time reaching some critical 



value that can be estimated utilizing Landau's criterium, T 



that de- 



termines a critical density when hydrodynamic evolution breaks up. Here, we 
will use a simplified version of this criterium and suppose the decoupling of 
hydrodynamical system when density in the center of the system reaches some 
critical value n f (typically, normal nuclear density) . Then the time t f , when 
the hydrodynamical evolution ends, can be estimated from the condition 

n(t f ,0)=n f , (27) 

If the hydrodynamical evolution stops before the condition (^6|) is violated, 
tf < t, then solutions (|2^) and ( p4[ ) describing quasi-one-dimensional expansion 
give complete hydrodynamical evolution, too. Let us find the conditions for 
such a situation. Supposing quasi-one-dimensional expansion and using (||) and 
l2|) we get 



Z(tf) 



no 

Uf' 



(28) 



Then using ( |24| ) we can find tf. If for t = tf the inequality ( pq ) is satisfied, 
then we can conclude that tf < t and the one dimensional expansion is valid 
approximation until the freeze-out time. 
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It is useful to give simple analytical estimations of the initial hydro dynamic 
conditions that guarantee quasi-one-dimcnsional expansion of the nuclear mat- 
ter. Let us suppose, for simplicity, that Zo= and hence u 2 = 2ji Z = Z , 
to = to in (|24|). Supposing that the upper time limit of quasi-one-dimensional 
expansion, t, is large enough so that 

Z 2 (t) > 4Z 2 , 

we can rewrite Eq.(|24|) in the form: 

Z 2 (t)^Z 2 {t)-3(Z Z 2 (t)) 2 / 3 . (29) 

After substitution (|2^) in ( p6| ) we obtain: 

3 %f-y ;< 1 v^^? , (30) 

I I i -y2 

It is easy to see from fl30P that for large enough ° zi " , quasi-onc-dimcnsionality 
can hold till the longitudinal scale becomes comparable to transversal ones: 



Z 2 (t) Z 2 (t) Z 2 (t) 



X 2 (t)+Y 2 (t) X 2 + Y 2 X 2 + Y 2 



cx 1 (31) 



It means that within time t < t solution (22), (EM) correctly describes the trans- 



formation of longitudinally compressed " ellipsoid" to " spheroid" form. Finally 
from d2g|) we get that Z(tf) < Z(t) oc ^/ Xfi + Yq and consequently tj < t if 



> < V3±R (32) 
rif Z 

Under such initial conditions one can expect that whole stage of the hydro- 
dynamical evolution can be correctly described by the approximate quasi-one- 
dimensional solution ( |22] ) and (pi[). 

Summary and conclusions: In this Letter we considered the time evolu- 
tion of fireball hydrodynamics describing an ideal gas, an elliptic initial density 
profile, a homogeneous temperature distribution and a Hubble-like flow distri- 
bution. For this case, the set of partial differential equations of non-relativistic 
hydrodynamics have been reduced to a set of ordinary, second order, non-linear 
differential equations, that can be solved efficiently by presently available nu- 
merical packages without the need of sophisticated programming. The initial 
conditions for these equations are associated with the initial elliptic sizes Xq, Yq 
and Zq, that could be linked with the overlapping geometrical sizes of colliding 
nuclei in heavy ion collisions and with the dynamics of the compression process 
during interactions during the pre-thermal time evolution. 

The general behavior of these hydrodynamical equations is determined an- 
alytically and related to the Hamiltonian motion of a particle in a repulsive, 
non-central potential. A first integral of motion has been found, corresponding 
to the conservation of energy in the Hamiltonian problem. It was utilized to 
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obtain an approximate solution for quasi onc-dimcnsional expansions and to 
determine the domain of applicability of this solution. 

The importance of the results is given by recent experimental findings in 
high energy heavy ion reacions, where various elliptic flow patters are observed, 
see refs. [[l0| |ll], [1^, |l3| for further details. In future studies, our results could be 
applied to gain insight into the interpretation of the above mentioned data and 
to describe nucleus-nucleus collisions with non-relativistic initial energies. Such 
conditions for a non-relativistic evolution may be reached in the mid-rapidity 
region near to the softest point of equation of state even in relativistic heavy 
ion collisions, if the pressure is not strong enough to build up a relativistic 
transverse flow. Due to the scale invariance of the hydrodynamical equations 
the solutions described here can also be utilized in other problems related to 
elliptic flows in non-spherical fireball hydrodynamics. 
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